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Abstract 



The Heisenberg antiferromagnetic spin chain with both dimerization and 
frustration is studied. The classical ground state has three phases (a Neel 
phase, a spiral phase and a colinear phase), around which a planar spin- wave 
analysis is performed. In each phase, we discuss a non-linear sigma model 
field theory describing the low energy excitations. A renormalization group 
analysis of the 5*0(3) matrix- valued field theory of the spiral phase leads to 
the conclusion that the theory becomes S'0(3) x 5*0 (3) and Lorentz invariant 
at long distances. This theory is analytically known to have a massive spin- 
1/2 excitation. We also show that Z2 solitons in the field theory lead to a 
double degeneracy in the spectrum for half-integer spins. 
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I. INTRODUCTION 



Antiferromagnets in low dimensions have continued to draw considerable 
interest in recent years, mainly because of their possible relevance to high Tc 
superconductors. A large variety of theoretical tools have become available 
to study them including non-linear sigma model (NLSM) field theories [1- 
6], Schwinger boson and other bosonic mean field theories fermionic 
mean field theories B, series expansions |TD[, exact diagonalization of small 



systems |TT[ , and the density matrix renormalization group (DMRG) method 



[|1^, |13l . In particular, one of the most useful tools has been the use of the 
non-linear sigma model (NLSM) field theories [1-6]. These theories display 
interesting features such as dynamical mass generation and existence of topo- 
logical soliton solutions, which are of relevance in spin models. In fact, it was 
the mapping of the 0(3) NLSM to the Heisenberg antiferromagnet which led 
Haldane to his celebrated conjecture that integer spin models would have 
a gap, contrary to the known solution for the spin-1/2 model. The exper- 
imental verification of this conjecture has given special impetus to the 
study of NLSM field theories for spin models. Ever since then, various gen- 
eralizations of the Heisenberg antiferromagnetic spin chain including models 
with dimerization and/or frustration |]T5| and coupled chains or Heisenberg 



ladders |]16l have been studied. 

In this paper, we study a general Heisenberg spin chain with both dimer- 
ization (an alternation 6 of the nearest neighbor (nn) couplings) and frus- 
tration (a next-nearest neighbor (nnn) coupling J2 ). Even in the classical 
spin 5* ^ 00 limit, the system has a rich ground state 'phase diagram', 



with three distinct phases, a Neel phase, a spiral phase and a colinear phase 
(defined below). For large, but finite S, a spin- wave analysis a la Villain [|18| 



can be performed which gives the spectrum of low energy excitations about 
the classical ground state. This is described in Sec. II. 

The Neel phase via the Haldane mapping to an 0(3) NLSM is well-known 
to be gapless for half-integer spins but gapped for integer spins [|I], 0. An 
interesting and important question to address is whether this distinction per- 
sists in the spiral and colinear phases as well. To answer this, we need to 
describe the long wavelength fluctuations about the classical ground state 
in these phases as well by non-linear field theories which are capable of sys- 
tematic improvement by renormalization group techniques. This has been 
identified in the spiral phase (for 6 = 0) @, |^ , and in the colinear phase for 
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5=1, which corresponds to two coupled spin chains [|I5], |T6[. However, while 
the Neel phase has been extensively studied, various aspects like the ground 
state degeneracy and the low energy spectrum are not well understood in the 
spiral and colinear phases. 

For completeness and to set the notation, we first review in Sec. Ill A, the 
NLSM and the RG in the Neel phase for arbitrary J2 and 6. For the spiral 
phase, in Sec. Ill B, we obtain the 5'0(3)-valued non-linear field theory and 
show using a (previously derived) one-loop /3-function that the field theory 
flows to an SO (3) x SO (3) symmetric and Lorentz invariant theory with an 



analytically known spectrum [|T9|. However, this does not directly lead to 
a distinction between integer and half-integer spins. But in Sec. HI C, we 
discuss how the presence of Z2 solitons affects the ground state degeneracy 
and the low energy spectrum and does provide a distinction. We are led to 
the conclusion that the system is gapped both for integer and non-integer 
spins, but the ground state is degenerate for half-integer spins (when there 
is no dimerization). In Sec. HI D, we briefly discuss the NLSM for a spin 
ladder. It is qualitatively similar to the one in the Neel phase. However, no 
topological term is induced and there is no distinction between the integer 
and half-integer spin cases. A shorter version of these results will appear 
elsewhere [0. Finally, we make some concluding remarks in Sec. IV. 



II. CLASSICAL PHASE DIAGRAM AND SPIN WAVE 
ANALYSIS 

The Hamiltonian for the frustrated and dimerized spin chain is given by 

H = J,Y. + i-^yS )Si ■ S.+l + J2 E ■ S,;+2 , (1) 

i i 

where Sf = S{S + l)/i^, the coupling constants Ji, J2 > and the dimer- 
ization parameter 5 lies between and 1. (We will henceforth set ^ = 1). 
Classically (for S — > 00), the ground state is a coplanar configuration of the 
spins with energy per spin equal to 



E, = S' 



y (1 + 5) cos ^1 + y (1-5) COS ^2 + J2 COs(0i + ^2 



(2) 



where 9i is the angle between the spins and S2i+i and 62 is the angle 
between the spins S2i and S2i-i. Minimization of the classical energy with 
respect to the 9i yields the following phases. 
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(i) Neel: This phase has 6i = 62 = n and is stable for 1 — 5^ > 4 J2/ Ji. 

(ii) Spirah Here, the angles 9i and 6*2 are given by 



cos 6*1 
and cos 62 



1 + 5 
1 

1-6 



4J2/J1 
4J2/J1 



+ 



4Jo 



1 + 52 

^ 4^ 
1-52 J, 



(3) 



where 7r/2 < < vr and < ^2 < ^i- This phase is stable for 1 — 6"^ < 
4J2/J1 < (1 - 6^)/6. 

(iii) Colinear: This phase (which needs both dimerization and frustration) is 
defined to have 61=11 and 62 = 0. It is stable for {1 — 6'^)/6 < AJ2I Ji- 
These phases along with the phase boundaries are depicted in Fig. 1. 

We now study the spin wave spectrum about the ground state. Before 
describing the details of the calculations, we state the main results. In the 
Neel phase, we find two zero modes with equal velocities. In the spiral phase, 
we have three modes, two with the same velocity describing out-of-plane 
fluctuations and one with a higher velocity describing in-plane fluctuations. 
In the colinear phase, we get two zero modes with equal velocities just as in 
the Neel phase. The three phases also differ in the behavior of the spin-spin 
correlation function S{q) = J2n{So ■ Sn) exp{—iqn) in the classical hmit. S{q) 
is peaked at g = {9i + 62) j^, i.e., at g = tt in the Neel phase, at ti/2 < q < n 
in the spiral phase and at g = 7r/2 in the colinear phase. Even for S = 1/2 
and 1, DMRG studies have seen this feature of S{q) in the Neel and spiral 
phases [|T3|. The colinear region has not yet been probed numerically. 

Since the classical ground state is coplanar (with the spins lying in, say, 
the X — y plane) in all the phases, it is more convenient to develop the 
spin wave theory by parametrizing the spin operators in terms of Villain's 
variables [|T^], rather than the usual Holstein-Primakoff variables [^. The 
spin variables at each site are expressed in terms of two conjugate variables 
Siz and (hi as 



St 

s- 

and Sf 



exp (20,) [{S+^y - (^Siz + ^y 



1/2 



1\2 



1\2 

S^z + 



1/2 



exp 



1 _d_ 

i dd), 



(4) 
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The periodic operator 0j (with 0j = 0j + 2tt) satisfies the commutation 
relation 

0i , ^ ] = ^ Sij ■ (5) 

In the large-S* hmit, the discreteness of Siz/S and the periodicity of 0j can 
be ignored. Si^ and 0^ can then be thought of as canonically conjugate and 
continuous momentum and position operators. In the classical limit — oo, 
Siz/S —>■ and is fixed at a value which is the angle made by the 
classical spin Sj relative to the x-axis in the plane. 

We now expand the square roots in Eq. and keep terms only up 
to order 1/5*. On expanding the angles 0j to quadratic order about 0j, we 
obtain 

Sr = exp m (S + ^) (1 + - \al) (l - ) , 

S- = exp (5 + ^) (1 - - (l - ) ' (6) 

where cxj describes in-plane fiuctuations (about the angle 0^) and Siz/S de- 
scribes out-of-plane fiuctuations. In the presence of dimerization 5, the unit 
cell contains 2 sites and has size 2a, where a is the lattice spacing. We 
therefore define 

Sn ^ = S2n+l,z 1 S^ = S2n,z ; 

and 0-^^) = CTsn+i , a^^^ = cxsn • (7) 
and their Fourier components 



and ai^^ = ^ ^ ai") exp(-z2A;na) , (8) 
where a = 1, 2, and 



[ al;^ , ] = t5'^'5,+,,,o (9) 
for the doubled unit cell. Substituting (||) in terms of the Fourier components 
of the doubled unit cell in (|l]), we find 

H = K + Sl'^^M'S^l + S' at^Bfa'^, ] , (10) 



k a,b 
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where < A: < 7r/2a and a,b = 1,2; and ^4^ and Bj^ are 2x2 hermitian 
matrices of the form 

_ f ai a2 + ase^^^" \ 

[a2 + a3e-^2^« Oi ) ' 

with oi = - Ji [ (1 + 5) cos6'i + (1 - (5) cos6'2 ] 

- 2J2Cos(6'i + 612) + 2J2Cos2ka , 

a2 = Ji (1 + 5) , 03 = Ji{l-5) , (11) 



and 

B 



\ b2 + bse-''^'"' bi 



with 61 = - Ji [ (1 + 5) cos6'i + (1 - (5) 0086*2 ] 

- 2 J2 cos(6'i + 62) + 2J2 cos(6'i + 6*2) cos 2A;a , 
62 = Ji (l + (5)cos^i, 63 = Ji(l - 5)cos^2 • (12) 

These formulae hold for all three phases; the angles 6i were given earlier for 
each phase. To find the spin wave spectrum u!{k), we have to diagonalize 
Aj^ and Bj., i.e., we have to find a matrix U (not necessarily unitary) such 
that det U — 1, and UAkW and U~^^ BjJJ'^ are real diagonal matrices 
with entries diag (ai,a2) and diag P2) respectively. Then the energies 
of modes of momentum k are given by = S^/al(3l and tU2k = 5'V«2/92- 
(There are two modes for each k because a unit cell contains two sites). 

In the Neel phase, both the energies vanish at A; = 0. The spin wave 
velocity c is defined to be the derivative du/dk at the gapless point. We find 
that c = 2JiSa^l — 5'^ — 4J2/J1 for both branches. In the spiral phase, we 
find one gapless mode at A; = (with a slope Co) and two gapless modes at 
2ka = 271 — $1 — 62- The latter are counted as two modes because the energy 
vanishes as one approaches that value of k both from the left and from the 
right; and the two slopes Ci are equal. The expressions for these velocities 
are quite complicated for general 5. (However, we find that the spin wave 
velocity Cg is always greater than ci). For 5 — 0, the expressions simplify to 

give Co = JiSa{l + U2/Ji)^Jl - jf/lQJl and Ci = c„y^[l + (1 - Ji/2J2)2]/2. 
Finally, in the colinear phase, one of the energies vanishes at = and the 
other at 2ka = vr. The spin wave velocities are equal at the two points and 
are given by c = 2JiSaJ {AJ25 / Ji + 5^ _ l)(2J2/Ji + 5 + l)/(25). One can 
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also check that in all the cases, the spin wave velocities vanish as we approach 
the boundary separating any two phases. 

The above calculation of the spin wave velocities c assumes long range 
order. For finite values of S, the spin chain models actually have no long range 
order, i.e., the two-spin correlation function goes to zero at large separation 
either algebraically or exponentially. However, it is still useful to calculate 
c because it is one of the parameters which enters the NLSM field theories; 
one therefore has a consistency check on the derivation of the field theories. 



III. NLSM FIELD THEORIES 



A. Neel Phase 



To study non-perturbative aspects, the NLSM approach is convenient 
since the RG can be used to improve naive perturbation results. The NLSM 
is well-known in the Neel phase and can be most easily derived as 

follows. The field variable is a unit vector 0. Since the classical ground state 
has a periodicity of two sites, we group the spins in two's and define 

T* S2i + S2i+1 



'2i+l/2 



2a 



and 02i+i/2 = 2^^'^^ ' ^^^^ 

These variables satisfy the identites 

f - = , 

1 a^^^ 

and 02 = 1 + ^ - ^ , (14) 

so that 0^ = 1 in the large- and long wavelength limit. (We will see below 
that i is proportional to a single space-time derivative of 0; hence long wave- 
length means that \ai\ << 1). We define the continuum space coordinate 
X = 2ia so that 52i,2i/2a = 6{x — y). The variables in Eq. (|^) therefore 
satisfy the commutation relations 

[ ia{x) , £[^{y) ] = i ea^y iy{x) 6{x - y) , 
[ 4(x) , (t)p{y) ] = i eai3-y 07(x) 5{x - y) , 
and [0„(x) , <Pp{y)] = 0. (15) 
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where a, /3, 7 = x, y, z, and 7 is summed over on the right hand sides. Thus 
i{x) can be identified as the angular momentum of the field 0, and £ = x vf, 
where vf is the momentum canonically conjugate to (p. The Hamiltonian 
(|l]) can be rewritten in terms of i and (p, and Taylor expanded in space 
derivatives, e.g., 

02i+5/2 = 02i+l/2 + 2a 02i+l/2 + 02.+1/2 + - • (16) 

For i, since it is already proportional to 0, we only need to keep terms up to 

^, /and where the dot and prime denote time and space derivatives re- 
spectively. Finally, we replace the sum J22i = J dx/2a to obtain a continuum 
Hamiltonian 



H 



I dx['-^{U^il-5)$f + ^J-], (17) 



where c is the spin wave velocity obtained earlier in the Neel phase, and the 



coupling constant = 2/ {S^Jl — 5^ — 4J2/ Ji). Thus large-S" corresponds to 
weak coupling. 

The Hamiltonian in (|T^ follows from the 'Lorentz-invariant' Lagrangian 
density 

• 2 

d> c0'2 ^ 



Here c = 2JiaSyl — 5^ — AJ2/ Ji is the spin wave velocity (a is the lattice 
spacing) and g^ is the coupling constant. The third term in (4) is a topological 
term with 6 = 27rS'(l — S). It is known that this field theory is gapless for 
6 = 7T mod 2tt with the correlation function falling off with the power 1 at 
large separations, and is gapped otherwise [0. For the gapped theory, the 
correlations decay exponentially with correlation length (, where ( is found 
from a one-loop RG calculation to he (/a = exp{27i/g^). Hence ln(^/a) = 



vrS'yl — 6"^ — AJ2/ Ji. For completeness, this is plotted in Fig. 1 for 5 = 
and 4J2/^i < 1- 
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B. Spiral Phase 



Recently, the NLSM for the spiral phase has been studied for 5 = |]. 
Since the classical ground state is generally not periodic, we will follow the 
treatment of Ref. 0. The ground state has 6i = 62 = = cos~^(— J1/4J2). 
The field variable describing fluctuations about the classical ground state 
is an S'0(3) matrix R{x,t) related to the spin variable at the i^^ site as 
(Si)a = SR^fjiib , where a,b = 1, 2, 3 are the components along the x, y and 
z axis, and n is a unit vector given by 

X cos z6' + y sin z6' + a£ . 
= — — r— — ^ y: ■ (19) 

I X cos ^c7 + y sm lU + at \ 

The unit vector fij describes the orientation of the i*^ spin in the classical 
ground state, and ai represents the small deviation from the classical con- 
figuration (it is not the angular momentum as in our discussion of the Neel 
phase in (^) ). The Hamiltonian in Eq. (|l|) can be expanded in terms of R 
and i, and Taylor expanded up to second order in space-time derivatives to 
obtain a continuum Hamiltonian 

H = J2 ^aMadb + tr {K'^gP) . (20) 

a,b 

Here M and P are diagonal matrices with 



M = diag (Mii,M22,M; 



33 J 



with M33 = 2J2aS\l + Ji/AJ2f, 
and Mn = M22 = M33[l + (1 - Ji/2J2)2]/4, (21) 

and 

P = diag {Pii, P22, P33) 



with Pu = P22 = J2aS'^{l - JiVieJf 



27' 



and P33 = 0. (22) 

The Lagrangian is then 

C = Ck-H, with Ck = aS{Te) -Y (23) 
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where 



T = diag (1/2,1/2,1) and Va = -l/2tabc{R^ R) 



be- 



(24) 



Ck is obtained by a path integral calculation using spin coherent states [pi]] . 
Finally, we integrate out the field C. which appears quadratically in C The 
resultant Lagrangian density is found to have an SO{?>)l x S0{2)ji symmetry 
and can be parametrized as 



^ tridtg'dtR Po) - I tiid.R'd^R Pi) 



(25) 



where c = JiSa{l + 4J2/ Ji] 
trices with entries given by 



1 — J1/I6J2, and Pq and Pi are diagonal ma- 



and Pi 



•^i^g 7T~2' 7T~2' ~2 ~ 7T~2 



diae 



The couplings gi are found to be 



9l 
9l 



9l - 

2gl 



and g\ = gl[l 



1 


1 1 


1 




2^' 7i~ 


"2^ 


1 


1 1 


1 




2^' 


"2^ 




'4J2 + Jl 






4J2- Jl' 




+ (1 


- J1/2J2) 





(26) 



(27) 



Perturbatively, there are three modes, one gapless mode with the velocity 
C92I g^ and two gapless modes with the velocity cgi/g^. Note that the theory 
is not Lorentz invariant because gig^ 7^ g2g3- However, the theory is sym- 
metric under 5*0(3)^ x S0{2)ji where the 5*0(3)^ rotations mix the rows of 
the matrix R and the S0{2)r rotations mix the first two columns. (To have 
the full SO{?i)l X S'0(3)i? symmetry, we need gi = g2 and g^ = g^, i.e., both 
Po and Pi proportional to the identity matrix.) The S0{3)l is the manifes- 
tation in the continuum theory of the spin symmetry of the original lattice 
model. The S0{2)r arises in the field theory because the ground state is 
planar, and the two out-of-plane modes are identical and can mix under an 
5*0(2) rotation. The Lagrangian is also symmetric under the discrete sym- 
metry parity which transforms R{x) — > R{—x)P with P being the diagonal 
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matrix (-1,1,-1). An important point to note is that there is no topological 
term present here (unlike the NLSM in the Neel phase) and hence, no ap- 
parent distinction between integer and half-integer spins. There is, however, 
a distinction due to solitons, as we will show later. 

At distances of the order of the lattice spacing a, the values of the cou- 
plings are given in Eq. (pT]) . At larger distance scales /, the effective couplings 
gi{l) evolve according to the /5-functions l3{gi) = dgi/dy where y = ln(Z/a). 
We have computed the one-loop /^-functions using the background field for- 
malism (Note that since the theory is not Lorentz-invariant, geometric 
methods cannot be used to obtain the /3-functions |Q.) The /3-functions are 
given by 



^{92) 
/3(53) 

and (3{gA) 



9l_ 
Stt 

9± 
Stt 

9j_ 
Stt 

9l_ 
Stt 



9l939i 



9l 



9i9i + 9293 



+ 25-15-3 



1 1 



~9l 



9l 



9M-2 
9i 

9l9i92 



- + 

92 



9l 



ably 2 

93 



9i9i 
2 1 



9293 



^9i93 ( — 

92 

+ 25153 



9l' 



9l 



9l 



+ 

„2> ^ 
54 



451^3 ( — 

9a 



(28) 



We numerically investigate the flow of these couplings using the initial values 
5j(a) given in Eq. (|27|) . We find that the couplings flow such that 51/52 and 
53/54 approach 1, i.e., the theory flows towards SO{?>)l x SO{?>)r and Lorentz 
invariance. Finally, at some length scale the couplings blow up indicating 
that the system has become disordered. At one-loop, C depends on J2/J1 
but S can be scaled out. In Fig. 2, we show the numerical results for ln(C/a) 
versus J2/ Ji for AJ2/ Ji > 1. Note that as ^J2/ Ji 1 from either side (the 
Neel phase for integer spin or the spiral phase for any spin), ln(C/a) 0, i.e., 
the correlation length goes through a minimum. Since 4J2/J1 = 1 separates 
the Neel and spiral phases, we may call it a disorder point. (For general 5, we 
have a disorder line 4J2/J1 + 5'^ = 1 and the correlation length is minimum 
on the disorder hne separating the two gapped phases.) 

The spiral phase is therefore disordered for any spin S with a length scale 
Since the theory flows to the principal chiral model with S0{?>) l x 5*0(3)^ 
invariance at long distances, we can read off its spectrum from the exact 
solution given in Ref. |]19[. The low energy spectrum consists of a massive 
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doublet that transforms according to the spin-1/2 representation of SU{2). It 
would be interesting to verify this by numerical studies of the model. DMRG 
studies |T^, |T3[ of spin-1/2 and spin-1 chains have not seen these elementary 



excitations so far. It is likely that these excitations are created in pairs and a 
naive computation of the energy gap would only give the mass of a pair. To 
see them as individual excitations, it would be necessary to compute the wave 
function of an excited state and explicitly compute the local spin density as 



was done in Ref. to study a one magnon state in the Neel phase. 

C. Solitons In The Spiral Phase 

Since the field theory is based on an S'0(3)-valued field R {x,t) and 
7ri(S'0(3)) = Z2, it allows Z2 solitons. The classical field configurations 
come in two distinct classes with soliton number equal to zero or one. If 
R(){x,t) is a zero soliton configuration, then a one soliton configuration is 
obtained as 

cos^(a;) sm6{x) 
Ri{x,t) = \ -sin6'(a;) cos6'(x) \ Ro{x,t) , (29) 
1 

where d{x) goes from to 27r as x goes from —00 to +00. (For convenience, 
we choose 6{x) = 27c — 6{—x), i.e., the twist is parity symmetric about the 
origin.) In terms of spins, this corresponds to progressively rotating the spins 
so that the spins at the right end of the chain are rotated by 27r with respect to 
spins at the left end. Since the derivative d^O can be made vanishingly small, 
the difference in the energies of the configurations Ra^x^t) and Ri{x,t) can 
be made arbitrarily small, and one might expect to see a double degeneracy 
in the spectrum. 

However, this classical continuum argument needs to be examined care- 
fully in the context of a quantum lattice model. Firstly, do Rq (x, t) and 
Ri{x,t) actually correspond to orthogonal quantum states? For the spin 
model, if the region of rotation is spread out over an odd number of sites, 
i.e., if the rotation operator is 

1 7r 

U = exp(^^ E (2^ + 2m + 1)5^), (30) 



n= 
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then RQ{x,t) and Ri{x,t) have opposite parities because under parity, 

m 

S't -S'_i and U Uexp{i2n ^ S^). (31) 



n=—m 



Since the sum contains an odd number of spins, the term muhiplying U is 
— 1 for half- integer spin and 1 for integer spin. Thus for half- integer spin, 
Ro{x, t) and Ri{x, t) are orthogonal and the argument for double degeneracy 
of the spectrum is valid. This is just a restatement of the Lieb-Schultz-Mattis 
theorem [^. For integer spin, Rq (x, t) and Ri{x,t) have the same parity, 
and no conclusion can be drawn regarding the degeneracy of the spectrum. 

An alternative argument leading to a similar conclusion can be made fol- 
lowing Haldane . We consider a tunneling process between a zero soliton 
configuration Rf){x,t) and a one soliton configuration Ri{x,t). (We choose 
coplanar configurations for convenience). Such a tunneling process is not 
allowed in the continuum theory (which is why the solitons are topologically 
stable) because the configurations have to be smooth at all space-time points. 
But in the lattice theory, discontinuities at the level of the lattice spacing 
are allowed. In terms of spins, this tunneling can be brought about by turn- 
ing each spin S-°^ in configuration Ro{x,t) to the spin in configuration 
Ri{x,t) by either a clockwise or an anticlockwise rotation. Assuming that 
the magnitude of the amplitude for the tunneling is the same (as we will show 
below), the contribution of the two paths either add or cancel depending on 
whether the spin is integral or half-integral. This is easily seen through a 



Berry phase |21| calculation. The difference in the Berry phase of the two 
paths from sf'^ to Sj-^"* is 27rS. Since the soliton involves an odd number of 
spins, the total Berry phase difference is mod 27i if S is an integer and n 
mod 2tt if S is half-integer. 

Now we have to check that the magnitudes of the amplitudes for tunneling 
are the same in both the cases. To see this, consider the pair of spins Sf 
and S^] which need to be rotated to Sj^^ and S^]. Since 6{x) = 2n — 6{—x), 
the magnitude of the amplitude for the clockwise rotation of S^^"* to S^^^^ is 
matched by the magnitude of the amplitude for the anticlockwise rotation 
of S^] to S^]. Hence, for the pair of spins taken together, the magnitude of 
the amplitude for tunneling is the same for the clockwise and anticlockwise 
rotations. 

Thus, tunneling between soliton sectors is possible for integer S (thereby 
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breaking the classical degeneracy and leading to a unique quantum ground 
state), but not for half-integer S (due to cancellations between pairs of paths). 
This agrees with the earlier Lieb-Schultz-Mattis argument. 

Although the NLSM model for the spiral phase was explicitly derived only 
for 5 = 0, we expect the same qualitative features to persist when 5 7^ 0, 
because the spin wave analysis shows that the classical ground state continues 
to be coplanar and there continue to be three zero modes (two with identical 
velocities and the third with a higher velocity. Hence we expect similar 
RG flows and a similar spectrum. However, the argument for the double 
degeneracy of the ground state for half-integer spins depends on parity being 
a good quantum number. When ^ 7^ 0, parity no longer commutes with the 
Hamiltonian and the argument breaks down. This is in agreement with the 
DMRG studies [|13[ (for periodic chains) which show a unique ground state, 
both for integer and half-integer spins, for 6^0. For open chains, the ground 
state is sometimes degenerate due to end degrees of freedom. To incorporate 
such effects, one would have to study NLSM theories on open chains which 
is beyond the scope of this work. 



D. Colinear Phase 



Finally, we examine small fluctuations in the colinear phase. The naive 
expectation is that the field theory would be an 0(3) NLSM, analogous to 
the Neel phase, since the classical ground state is colinear. We can show 
this explicitly for 6 = 1 which is called the Heisenberg ladder [^]. The 
field theory in this limit can be derived using the classical periodicity under 
translation by four lattice sites, similar to the derivation given in Se. HIA 
for the Neel phase. For a set of four neighbouring spins, we define 



(x — a) = , i {x 



2S ' ' ' 2a 



and {x + a) = ^^-+3-^84^+2 ^ ^ = ^^-+3 + 84^+2 ^^^^ 

where x = (4i + 3/2)a is the mid-point of the set of four spins. We then write 
the Hamiltonian in terms of the fields (p and £, and Taylor expand to second 
order in space-time derivatives to obtain the Lagrangian (|1^) without a topo- 
logical term. We now find c = AaSJ J2{J2 + Ji) and g'^ = ^ J {J2 + Ji)/ J2- 
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The absence of the topological term means that there is no difference be- 
tween integer and half-integer spins and a gap exists in both cases. In fact, 
the NLSM predicts a gap for any finite inter-chain coupling, however small. 



This is in agreement with numerical work on coupled spin chains ||26 

IV. DISCUSSION 

In conclusion, we emphasize that this is the first systematic field theo- 
retic treatment of the general Ji — J2 — S model on a chain. Although all 
experimental spin chain systems known to date, like NENP and Sr2Cu03 , 
are in the Neel phase p7[1 , it would be interesting to find an experimental 



system with sufficient frustration and dimerization to probe the spiral and 
colinear phases. These phases could also be studied using numerical tech- 
niques like DMRG. The field theoretic treatment of the spiral phase leads 
to the interesting possibility that the low energy excitations of integer spin 
models may be massive spin-1/2 objects. This again is a possibility which 
could be looked for experimentally or verified by numerical simulations. 

The field theories for general 6 in both the spiral and colinear phases are 
still open questions. Although the results are qualitatively expected to be 
similar to the 6 = case in the spiral phase and the 6 = 1 case in the colinear 
case, quantitative features such as the dependence of the gap on the coupling 
strengths will require the explicit form of the field theory. 
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Figure Captions 



1. Classical phase diagram of the Ji — J2 — S spin chain. 

2. Plot of \n{C/a)/S versus J2/J1 for 6 = 0. For AJ2/J1 < 1, ln(C/a) is given 
by the one-loop RG of the 0(3) NLSM for integer spin. For 4J2/J1 > 1, 
ln(C/a) is given by the one-loop RG of the S0{3)l x S0{2)r NLSM. 
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Fig. 2 
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